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I XMITODUCTIOM 

In our tvjo previous reports attention was paid to tlie compilation of 
cloud photointerpretation keys for Landsat imagery so tlnat the recognition 
of various cloud types might be pxit on a firm footing (Barrett and Grant, 
1976a) , and eyeball and macliine-assisted coirparisons were made betvreen 
estimates of cloud cover imaged ty Landsat and observed from the ground 
(Barrett and Grant, 1976b). This latter work included the development 
of an appropriate methodology, and the compilation of sone initial results. 
Jn due course it is intended that these results should be suppl^aiientwi by 
others for that partial of the study period for whion imagery had not 
become available (generally February and March, 1976). Since we had not 
received our flrial consi,gnments of images by mid-August 1976 the ccaipletion 
of that part of oijr exercise has still to be effected. Table 1 and 
Figure 1 indicate the frames received since the preparation of the Third 
Quarterly Report: Figure 2 displays the distribution through time ■of ou).-> 

entire archive of Images as it stands at tlie time of writing'' 

It was thought convenient, whdJLst we were ai‘ 7 aiti'*S the remainder of 
the images far our study period, to experiment with a furidier section of 
our work which seeJes to evaluate Landsat aloud images in coHparison ivdth 
those from other satellite families. This section is concerned with 
simultaneous analyses of data in all four Landsat multispectral sensor 
wavebands, leading to an appreciation of the suitability of tlie Landsat MSS 
wavebands for clowrl analysis and idcntifijcation. Many studies liave con~ 
sidered tlie multispectral responses of Landsat-vxewed surface-features; 
many more have considered the significance of multispectral data for 
differentiating both between and within classes of a wide variety of 
surface features (see e.g. NASA, 1975). However, we are unaware 
of detailed studies of the multispectral characteristics of Landsat-viewed 
clouds, tdiough Danko (1974) addressed himself In par t to this question. 
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LANDSAT COVERAGE OP THE BRITISH ISLES; 

Tabulation of Individual Frames 
(see also Table 2 in Barrett and Grant, ^975(b), Table 1 in 
Barrett and Grant 1976(a), and Table 1 in Barrett and Grant 
1976 (b)). 
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Fig* 1(b): Landsat coverage of the British Isles, Cycle 18 

17th December, 1975 - 3rd January, 1976* 























Pig. 1(f): Landsat coverage of the British Isles, Cycle 22 

27th February - 15th March, 1976, 
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Indeed, there seems to have been a vjidely- -accepted belief that Landsat 
ijiBgery is ill-suited to cloud studies because cloudiness is thought to 
saturate tlie pictures easily (see, e.g. Danko, ibid. ). 

Dependent as we have been on image dala, not CCTS;. we are not in a 
position to comment conclusively can this belief: there are more 
uncertainties in the analysis of Landsat images -tdian CCTs, stemming mostly 
from picture processing. Taking as many of these uncertainties into 
account as we are readily able to do, and working at something less than 
tile ncminal full resolution of ttie Landsat inages, it appears that the 
principal families of clouds do possess distinctive spectral signatures, 
and yield different histograms of picture brightness over mesoscale areas 
of cloud fields. 

Tlie inrplications of these findings are discussed in relation to 
otlier work in this Department concerned with automatic (objective) 
analyses of cloud imagery fror. meteorological satellites. The biggest 
potential advantage that an automatic multispectral cloud mapping 
mstliodology might enjoy over one based on single waveband data seems to be 
that that might be able to utilise tiie full spatial resolution of the 
original observations, rather than the reduced resolution which results 
from most methods of texture analysis. 

The choice of sampling ar-ea was important in view of the characteristics 
of the cloud fields we wished to investigate, the imagery we chose to use 
in our analyses, and the capabilities of t3ie equipment available for the 
task: 

a) Cloud field characteristics. Althoii^i it .ould have been 

possible to select a brightness threshold to represent the edge 
of the cloud portrayed in any image so that non-cloud areas 
might have been excluded from the study, we did not think it 
necessary or indeed desirable to take such a step. We did not 
w’ant to make initial assumptions which might influence the 
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shapes of image brightness/frequency distributions, especially 
since tlie characteristics of cloud margins, and 1die spaces 
between adjacent structural elements of cloud fields mi^t 
significantly assist in some aspects of the autonetic recognition 
of different faniilies of clouds. 

.,) Tlie form ox the imagery. In viev7 of the added problems that 
would liave been associated with copying and enlarging the 70 mm 
transparencies provided by I'lASA, we decided to analyse these 
images themselves. At this scale 1 min represents approximately 
3.37 km on the ground. 

c) Tl^ie microdensitometer . The instrument we used vjas a System P-IQOO 
Photoscan , manufactured by Optronics International Inc. This is 
a rd.^“speed digital microdensitometer. The machine incorporates 
an electro-optical x'otating drum which converts photometric data 
on film negatives or positives to digital form for computer 
processing. In order tlsat an image may be scanned it is placed 
ovear an opening in the drum, and claiiped to it so that the film 
adheres exactly to its madiined cylindrical surface dimensions. 

In the present study a sanpling aix^a of 10 mm^ was chosen. VJe vdshed 
to sanple a sufficiently large area to establish whether the brightness 
characteristics of fields of different categories of cloud are indeed as 
clear as eyeball investigations seem to suggest. It seemed likely that 
cloud xi.elds should be distinguishable from one another if viewed 
sufficiently broadly to account for their textural cliaracteristics j 
whether cloud faniilies could be differentiated on the basis of their 
reflectance characteristics per unit area seemed to be mere doubtful. Our 
hope, was that, by compiling sufficiently large populations of summary 
statistics for quite large areas, both issues might be elucidated simulta- 
neously, the first throu^ areal brightness/ frequency distributions., the 
second through the spectral reflectance graphs constructed for each 
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cloud category from the results for the modal classes in each briefness/ 
frequency histogram. 

Accurate registration was attenpted for each sot of ii.Tagery using a 
clean perspex mask on which were placed tlic locations of the four regis- 
tration marks of each frame. In tliis way -1000 four images of each cloud 
scene coiild be re;gistared and a common san 5 >ling area identified, 

II TCaMIQUES 
1. Selection of imagery 

It was established earlier in this project that major cloud types can 
be differentiated and identified wi1h confidence in Landsat inegery. The 
cloud families recognised thus were cunTulordmbiforrn. cumuliform, stratiform, 
stratocumulifom ernd cirriform. Certain Landsat image: were selected to 
represent "classic" forms of tliese cloud families and tliair chief members 
or cloud types . (See. Barrett and Grant, 1376(a)., especially Mates 1-5). 

For our multispcctr'al study we cliose one example of each cloud family, 
plus an cxai.ple of tiic cumulocongestus cloud type, the latter on tiie 
grounds tbat cumulus humilis or irsdiocris and cumulccongcstus have sub- 
stantially different appearances in the Landsat Images. Since all four 
spectral bands wei>! to be investigated, there was a set of 24 cloud images 
to be analysed. 

All 24 iiTHgGs depicted sea areas around the British Isles. Background 
brightTiCss constitiites a serious problem for any programme of automatic 
analysis of cloud ii.nagery. Methods have been devised to solve this 
problem (see, e.g. Hiller, 1971), but none have been (or could 
be?) entirely satisfactory. Vie Mslied to reduce, rather 'ttian solve, the 
problem by analysing clouds viewed-over surfaces with relatively even and 
constant bri^tness responses. Sea surfaces meet these conditions much 
better than land surfaces, especially in the Britisli region, where con- 
vergences of major currents do not occur. 
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2. Selection of ganpling area 

This was defined on tf.e mash by scraig>it -edged, irachine-cut, metal 
foil strips which overlapped to provide a precisely measured imaging area. 
In operation each image was orientated correctly on the nask provided for 
its cloud category, oj'id tlie image vdth its mask wexx; then located on the 
imaging drum. 

The Koehler’ illuinination system of the Fhctoscan ensures uniform 
illununation and focus sir^g on tlie film surface via • turret mounted apertures. 
Tiie lig.ht trarismi'tted through eacii image is measured using a photo- 
detector > ai'id converted to 256 grey-levels . 

The choice of scaijrdng apei^turvc is important . as this determines the 
sampling size for the irragery. Ttie aperture siz<cs available for imaging 
on the Bristol Urdversity machine are threefold, namely 25, 50 or 100 ym 
square. In the R‘\SA Landsat Data Users Handbook (1B71) (see H 1.4, 
''Photographic Hicrc-Iiaage Quality’ ) the folicwing is found; 

"Because of MTi’, granularity, and sensor and recording systems 
considerations, scans using aperturiac sxiallcr than - 20 pm 
diameter v.ill be essentially meanir,gle3s. scans witii aperture 

sizes of - 40 piu v.ill probably not correlate wdl with macro-density 
readings even ntcri assun’iin>7, tiiat the inwistigator has made necessary 
corrections fremi specular to diffuse densit;/'' 

Vs/ith such considerations in mdnid, tlic largest iatiging aperture size 
available (100 pm) was selected in order to miimise the sources of error. 

An ^>erture of 100 pm square corresponds on the irrage to an area, of 
approximately 337 iriotrso square on the ground. Thus approximately 16 pixels 
were analysed in each of our microdensitaT.eter spot readings. 

3. Operating practices 

The illuminating and imaging optics ox t}ie Photoscan are mounted on 
opposite sides of a carriage tiirough v;hicli the cya'lnder drum rotates. 


21 , 


The optical density of eadi image is measured every 100 yra along the 
circumference of the drurii (Y direction) vjithiri the pt^-detemirted area. 
After each revolution the 'C* carriage is stepped in the axial (X) 
direction by 100 ym. In practice these processes v;ere repeated until 
the whole 10 mm^ area of cad'i image had been scanned. Once per revolution, 
through an opening iii the drum opposite the filn position, ttie densitemeter 
photo- -detector system is reset to a given Jo'iOwn value vJldch represents 
im optical density of zero as defined by the air path tlirougli tl-ie slot. 
Since the drum sijeed is high (8 revs, per second) the drift of the 
instrument is far less significant th»ci'. the Iccist significant bit of 
derisitjf data. 

Hie dotootor voltage resulting from light tx'aiisrnitted tlirough each 
image is anpiificd logarithmically, digitized, displayed, and recorded 
by interfacing the scanricr v?ith a magnetic tape recorder. Each number is 
representative of a grey-lcvcl in thic selected density range. In our (sase 
the density range was 0 to 2D. VJithixi t]iis range density bears a linear 
relationship vn.th the outpUT gray levels. 

Optical density (D) is defined as: 

D = log ^ 

x-jhare 

li is the light intensity’ inpinging on tlie. detector through an 
air path and 

It is the li^t intensity of the tr<3nsmitted light. 

In addition to tlie 10 mnr sanpling ar-a on each image, the 15 step 
grey-scale tablet vas scanned. Ihis scale lias ledoc'gone the same copying 
and prcjcessing as the image to which it is attached (NASA, 1971, pages 3-5) . 

4. Data output 

Output from the microdensitometer was recorded on a magnetic tape 
unit for subsequent prvTcv^ssing on a PDF 11/45 mini-c^’:ptrter. 
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Initial processing vras concerned with tlie recording brightness 
values of the 15-step grey-scale tablet of eacli image. A number of data 
records were exajmLned for each image and the upp^sr boundary value of each 
of the 15 steps was determined by inspection ox -tlie ccrfputer print-out. 
Preferably sane tyf-e of edge-detection logic (e.g, Roserrfeld, 1970) 
should have been eirployed., but tliis time did not permit. Once the upper 
boundaries of each stop had been deterrxiined these values were fed back 
into the computer in order that the frequnT-ici,' distribution of briefness 
values fer each image could be determined. This prcccdurc was adopted 
so tliat we might be able to relate different frequency classes to actual 
tones on each Landsat immge. rather* than some arbitrarily chosen density 
value which might be inappropriate in viev7 of the apparent variations in 
image processing. 

Besides the frequency distributions of the brightness values, other 
summary statistics were coiiputed also, including means, standaitl 
deviations, skewness and kurtosis. 

The mean brigJitiiess value for each imcige was deteri^nined in the foUo^^ 
ing manner. The store size of the conputer is limited (24 k) and was not 
large enough tc Jiold all the values of each image ( 8000-9000) . Thus the 
mean brightness value of cacli record (Y-ioas scan) was calculated and 
the final niean ms taken as tl;e mean of the record means. The Standard 
Deviation v;as oalculatfed using tiie final mean derived as outline above. 
Skewness and kurtosis values v;ere not found directly on the PDP 11/45 
because of problems with the computer program. These V7ere calculated from 



'the frequency distribtrtions using a Hewlett Packard 9801A calculator and 
an appropriate library program (Prog. l*-3). 
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III RESULTS 

1. Brightness Gharactferlstics of l,andsat positive transporcaicies 

Fcet analytical purposes , all the biughtness values produced by the 
microdensitometer in terrns of 256 grey- levels vjere converted to optical 
density vcilucs in tiie rarige 0 to 2D. Tliis was a straight-forward 
procedure-- as the density bears a linear relationship with the grcy-levelSj 
as noted previotisly. 

Fig. 3 shosfjs the 15-sTep grey-saile tablets, plotted as a function 
of density and waveband. Tlie divisions correspond to the upper lindt of 
each step, the procedure for detemdning wlTidi was outlined above. Step 8 
of each wedge is shaded to act CiS a reference. Step 1 (white) is at the 
left-hand side, and step 15 (black) is at the r'ight-hand side. 

Perhaps the most important fact emerging from this diagram is the 
wide variation of step vri.dtli. Generally, the bri^ter steps, with densities 
less than )D, ara narrov 7 er' than the darker steps with densities greater 
than LD. The first stop is somewhat anomalous in tiiat it has no lower 
limit. Another notable feature is the variation from image to image in 
position of similar steps. Hds variation occirs both bet'woen vjavebands 
for a particul:Sr frame, rend also betv/een one frarne and another. ThvB, for 
example, in the illustration of cumulus humilis, corresponding steps of 
hand 7 are generally triglTtor 'Idian those of tand 6. Sucli tUfferances are 
due in part to variations in the processing of Landsat images. 

A serious problem of interpretation therefore ai'ises. Are the 
differences in briglitness in a particular sceriC in different v/avebands 
due principally to differences in target iv?flectance, or differences in 
processing? It is clear from our results that tlie influences of image 
processing must be borne in mind constantly as we move towards dravjing 
conclusions from tlicse data. Ideally, the tt-jo conponents of variation 
should be deteradned and separated i this vje could not adiicve in our 
study. 
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Fig. 3: Grey-Scale Tablets. 


25 


2. Bpjfghtnssa diaracteristxos of six categories of clouds; smphical results 

Figs. 4 a to f sliov? frequency distributions of individual cloud 
types in each waveband in relation to the 15 step gr\ 2 y“scale. Frequencies 
liave been converted from number coui'ts to percentages to facilitate com- 
parison between different. cloud types. Notwith';:^ ending our earlier 
ccmiiuBr.ts on variations in picture^ processing, some inpoortant differences 
do occur. 

a) Cumulonimbiform (Fig. 4(a))'. Bands 4 and 7 show apparently 
similar distrib’utions .j with band 7 being displaced sligiitly 
towards the darker end of the scale. Bands 5 and 6, however, 
possess a strong modal category, well displaced towards tiie 
brigfitei-' end of Id'n scale. 

b) Cumulxis congestus (Fig. 4(b)): The frequency distributions of 

all four wavebands range across the whole step-wedge scale. 

All the dis'tributions are essentially similar, seme minor 
differences, occurring at the dark end of the scale in steps 
14 and 15. 

c) Cumulus humilis (Fig. 4 (c) ) : Unlike the cumulus congestus 

graphs of Fig. 4(b) , tliese distributions are restricted in range 


to the darker end of the scale. Tiiis reflects -the lower bright- 
ness responses of tliese small cloud cells and the larger areas 
of dark sea background between t'nem. 

d) StratifaciTi (Fig. 4(d)): Each waveband sliows a narked clustering 

cf 'values at tl'.e bright end of the scale, indeed in band 5 
a3mo.st 100% of the data points are contained vn.thin the brightest 


step. Tliis is undoubtedly due to extensive saturation of tlie 
image by this cloud type. The distributions of bonds 6 and 7 


shov; tails extending t«x;ll into the dorker portion of the scale. 



Stratocuniuliform (Fig. 4(c)) ; 
tributionn, extending across 


Bands 4 and 5 have similar dis- 
wide range of values, peaking 
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Band 7 

N.B. Values <r/o not shown 

Pig- 4(d): Stratiform: Frequency Distribution, 




NB. Values <Vio not shown 

Fig* 4(e): Stratocumuliform: Frequency Distribution 
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gently in tiiS centre of the scale. In hand 6 there is a 
sligtkt displacGiEcnt towards the dark {rricl of the scale. In band 
7 lixere is a sigrdficant displacenient towards the darker end of 
the scale;, and the range of values is more restricted, 
f) Cirrdform (Fig. 4(f)); Bands 4, 5 and 6 are similar to one 

anotl'ser, x-jith distributions located in the darJeer portion of the 
scale ctnd a fairly wide srange of values (over 8/9 steps) . Band 
7 possesses a marked displacement tov;ards the dark end of tiie 
scale and the range of VetLues is limited to the darkest 6 steps. 
Figs. 4 (a) "(f) reveals therefore 5 that differences in brightness 
occur both betvjeen the x-javebands for individual cloud types ar vfell as , 
more strikingly, betv.’aen similar wavebands for different cloud types. To 
facilitate the examination of these differences ^ cumulative frequency 
distributions v?sre plotted. (See Figs. 5(a) -(d)). 

These graphs consist of plots of percentage cuiuulative frequency on 
an arithmetic probability scale , against image densdty on an aritlimotic 
scale. First, the average density of eacli step was found. In the case 
of step 1 the average density wejs found by detarudiiing tlie average density" 
range of steps 2 to 5 inclusive. This value vras then lialved, and sub- 
tracted from the upper limit value of step 1 to provide an average value 
for the bri^tness of that step. The advantage of using aritlimetic 
probability paper in this analysis is that if the distributions plotted 
are nonnal, then a straight line will occur on the graph. A brief exanina- 
•Idon of these diagrams shows tliat the distributions xvGre not noimal.. The 
chief points that emcarged from tliis analysis Xvwra as follows; 

a) Band 4 cumulative frequcnc/ distributioiia (fig, 5(a)). Here 
stratiform ard cmuloniMjifarm clouds appeal’ to give siiailar, : 
results, both having the majeritx' of pheir Vcilues at tLie bright 
end. of the . scale » Similarly cirriform and stratocurniilifom clouds 
are similar to each other, their curves lying close together 






acaross mudi of the graph. At first si^t this may appear 
strange as the graphs of •tiieiu:' respective frequency distributions 
(Figs. 4(g) and 4(f)) c^percntly ixvealtxl ccansiderable differences. 
However j one examines Fig. 3, it Ccin te seen that strato- 
eunulifom gsrey^saale values are considerably bri^ter liian 
GGfrrssponding cirriform va3.ues. This is taken into account ncre 
fully in Fig. 5(a) where actual density values are eiqf>loyed on 
tile abscissa* not grcy-seale step numbers as in Fig. 4. 

Generally, most cloud types follow paths wiiich are ropgiily parallel 
over much of the scale in Fig. 6. The cest notable exception to 
tliis is c\3nulocongestus . The p>lot for this cloud category starts 
wall to the left of the grapli and generally rises less stet^ly 
than those for the other cloud types exc^t at the dark end of 
the scale. Thie tendency for this cloud type then* is for its 
plot to cut across the plots for cirriform, stratoomulifcam 
and evanulus humilis. 

b) Baik! 5 cumulative frequency distributions (Fig. 5(b)): These 

are essentially similar to the band 4 plots in Fig. 5(a). 

However, the siaratiform curve is shorter due to tlie clustering 
of its values at the bright end of tdic scale. Tha eumalos 
Gcsngestus curve again lies across those for the other types. 

c) Bcffd 6 ciJimLative frequency distributLcns (Fig. 5 (c)>; For 
the first time stratiforTii ext^ds into dailcer pedrtiqns of -die 
stsde Ihsm dQe& a 

inerkedLy displaced tewards Ihe dar in 

^ " "6oS|jariSESi -itj^ ■ 'CumuluS^^^^ \ 


tcxrards darker suLthou;^ it reties its 

" j^siiel ^ ^sit^on' with: 'fespeerL'' to ' diiriisEin 
Cumulocongestus intersects the other distributions at quite 
I, as in bands 4 and 5. ths dafk;^^ 


' Ki^P:EOiJ^UeiiM.Lir i (;F;. ib 



d) 



I 

■ 38 .' ■ ■ ■ ■ i 

• ■ S 

i 

i 

the final 30% of its distribution follows a similar path to liiat I 

! 

of cumulus hunalis. \ 

Band 7 cumulative freqiBnGy distributions (F%. 5(d)): In this | 

waveband, cumulonimbifom has been displaced quite jftarkedly 
towards tlie dark portion of the scale. The wliole of its curve i 

j 

is on the darker side of the stratiform curve. CiiTriform and i 

stratocumuliform distributions now foHd<r very similar paths, I 

both having been displaced nirirkedly towards the darker portion of 
the scale than in band 6. Cumulus humilis is darker tiion in 
band 6. Tlds is probably due to the darkening of the sea back- 
ground which ecmprises a larger proportion of this inHge between 
the small Gumulus colls . This effect is also seen in the case 
of cirriform and stratocumul iform, both of vdtich contain sdne 
areas of seit in the imaging area. 


These diagrams , Figs. 5 (a) -(d), show tiiat marked differences in 
reflectance characteristics do occur betvjeon cloud types and that these 
diffcrenc'GS change frcjn ont; waveband to anotiicr. It is thought that such 
differences are indeed due principally to t,argct characteristics, not 
processing practices , for more account is taken of the grey-scale differences 
in these diagrams through; the use of the actual density values of the steps 
rather than tJieir numbers in the grey-scale. 


3. i^i^tness eharacteristios of six categories of clouds . suntnary statistical 
results 

In addition to such graphical descriptions of the bri^tness frequency 
distributicais, it is possible to gain further iiisight into the nature of 
cloud differesnees by deriving suntnary statistics via the msldiod of monents. 
This is a ccxiputatiGnal method; the analysis of the mean and standard 
deviation was carried out on the PDF n/4S nuni-conputer'. Skewness and 
kurtosis coefficients were determined on a Rewlctt-Packard caiaxlatbr 





fccin tJie frequency <3istributions derived by the PDP 11/45 . The results 
of the ccraputations tire presented in Table 2. 

•Jliere are some interesting differences between I3ie sunanary statistics 
and the inpressions gained of the distributions presented in Figs. 4(a)-(f) 
For example, the cujnulcniihbifonn statistics shovj tlTcVt the kurtosis 
("peakednesS'') of the distributions increases steadil y frem band 4 through 
to band 7, exandnation of Fig. 4(a) alone rc>uld probably have led one 
to conclude thiat bands S ^d 6 were more ’’peaked’’ than either bands 4 or 
7. Cumulocongestus clouds shov? relatively large standard deviations, a 
fact less mexpected when Fig. 4(b) is examined. A snail negative 
skevmess value is noted in band 7 which corresponds to the increased 
frequency of values at the dark end of the scale (step 15). The cumulus 
humilis distributions (Fig. 4(c)) are asymmetrical in nature, and; tipis 
produces negative skevzness values for all wavebands. Mean bri^tness 


values shov; a similar trend in all. the cloud types. Generally, the mean 
brightness of each successive waveband is slLj^i’tly darker than the 
preceeding viaveband. Tliorc are two exceptions to this general rule. 
j3otn occur in v;aveband 5 in the cirriform and stratifom cloud types. 


Such a decrease in brightness vri.1±; waveband is probably explained by 
tlie decreasing bri^itness of the background area (the sea) with successive 
wavebands . 

Besi'ies this slight ancsiialy in the general trend of mean brightness 
values , stratifciEnn clouds also present some other notevrorthy Sunmary 
statistics. The skewness values recorded are very large, indicating 
extreme asymmetry of the distribution. This, is perhaps not surprising 
Vivien Fig. 4(d) is exandned, especially with respect to band 5. However, 
the kurtosis values seem excessively high wh ?n one considers that in 
sediment analysis values of kurtosis greater than 3.(X) indicate extreme 
leptokurticity (peakedness) and are uncommon. 









Spectral rcflect.:jic£. grapiis for clx c:nt3i.xcries of clouds 

Tlie results siamzarised ii*i sections Cl) - (3) above are all related 
to characteristics of cloud fields viewed over areas of about 33.7 km^ 
on the .ground. It would be interesting to investigate the degrees of 
similarity between those results and others for a range of smaller ai'eas. 
Whilst we Eiust be cautious not to read toc> much into the results obtained^ 
it seems legitimate to attempt to inteiTpret them in terms of what they 
may indicate concerning tlie cdiaracteristics of cloud elements viewed by 
the microdensitometer over areas of about 337jii2 on the ground. Taking 
the modal classes in Figs. M-(a)-Cf) to indicate ti'ie most representative 
brightn^s vLiues for each cloud category ii eadi waveband we have compiled 
a set of spectral reflectance graphs for the six cloud examples (see Figs, 
6(a)^‘Cf)), It appears therefrom that the spectral reflectance signatures 
for different categories of clouds possess some marked differences both 
in form end position on the exes. Not surprisiigly the eumulifcm and 
cumulocongestus curves .3re generally similar, whilst the stratiform and 
coBraiLonimbiforTO curves are disappointingly alike . 

It would seep .1 w<3ll wrthwhile investigatil:ig the multispectral 
differences of Ljoidsat- imaged clouds more carefully j pvjrltaps througii the 
use of digital data from the CCTs so that the full resolution of the 
imagery might be esi^loyed ahd picture processing problems thereby 
eliminated^ 


Discussion 

It is in the field of sediiBsntacy petrology and geology tliat this lype 
of statistical analysis has reached its groatest sophistication. The tech- 
ruLque of moment sjoalj^is heis been used to ide^itify and differentiate sedi- 
ments (e.g. Folk and Word, 1957) and in the iiiterpretation of the 
deppsitional histories of different sodinisgts Ce.g. Greenwpe3d^ 1972).^ 
Generally j bivariate plots have been employed ^ with eadi of the four 
sumrrmry statistics being plotted in turn against each of tdie others. . , ^ 


Cumulonimbiform 


Cumulus-cpngestus 




Cumulus-humilis Stratiform 




StratoGumuliform Cirriform 



SR = Spectral reflectance 
D- Density 



Fig, 6: Spectral Ref lectaaice Graphs, 




43 


However, Folk and Ward (1957, p, 23) inanaged to produce a four-variate 
graph to represent the relationships between the four sunmary statistics 
in One diagram. Often, the analysis of such data has been extended by 
using clustering and discriminant analyses as classification procedures 
to allow new data values to be assigned correctly to previously determined 
groig>ings. 

It is felt that this latter type of analysis may prove most useful 
in attempts to assign new cloud information in Landsat data to previously 
determined cloud types, defined in terms of their momeint statistics. 
However, some problans would have to be solved before this were to be 
possible. For example, our results for stratifonn clouds indicate that 
the application of such techniques to non-normally distributed data without 
sane prior transformation can . lead to misleading, possibly erroneous, 
results and conclusions. One solution might be to transform the data 
after initial statistics had been calculated. For exanple, Folk and 
Ward (1957) produced a transfonned kurtosis statistic (K*) derived fron 
the ordinary kurtosis value, (K) . The relationship is as follows: 



.K 

"(K + 1) 


If this trvwsformation had been applied in our study, the range of 
kurtosis values wuld have ranged frcm about 0.5 to 1.0, much less than 
the range we found. This transformation produces an approximately 
. normal kurtosis distribution (Folk and Ward, 1957 , p. 15). 

A further problem with our distributions is that at the brightest 
end of the scale they ;^ "open-ended"; no lower, limit was defined. 

Folk (1965) , in discussing grain size analyses of sediments, felt that 
the application of the method of ncments to such open-ended distributions 
mi^t not be justified. This may be pertinent to our extreme results 
for stratiform clouds, where likely saturation of the image caused 
many values to be located in the unbounded step 1 of the grey-scale 
step-wedge. , ■ -- 


It is apparent, therefore, that considerable work rennins to be 
done if a wortl^ aim is to produce an operational cloud identification 

scheme for Landsat-type imagery based on cloud briglitness statistics. 
However, the difficulties should not prove insurmountable if careful 
thou^it is given to ejqjerimental design and Operation. 

ACCQMPLISHMEtirS 

Analysis of san^ile multispectral images for six different categories 
of clouds have revealed marked differences between the reflectance 
characteristics of clouds and cloud fields or cloudy areas in the visible/ 
near infrared region of the electromagnetic spectnm. Although the 
results obtained so far relate only to clouds imaged over coastal waters 
of the British Isles in summer the method by which 1he results were 
obtained could be applied quickly and cheaply to the analysis of further 
nnages. 

PROBLEMS 

Several problems of note were encountered in this stage of our 
EKjCS Follow-on Programme Study. These included: 

Processing problems 

These were of ti-jo different kinds, namely those already discussed 
at sane length related to evident variations in the processing of Landsat 
imagery by NASA, and those related to the equipment available to us for 
image analysis. The two are interconnected in the sense that the poor 
quality of the copy negatives provided by NASA (see Barrett and Grant, 
1975a P.8 ) determined that we shovild use the 70 mm positive 
tr’ansp^rencies obtained from the same source for the microdensitometry 
investigation, despite Iheir anaH size for this purpose, rather than 
enlarged positive transparencies generated from our archive of 70 im 
negatives. Since the Optronics Fhotoscan has a limited nunber of spot 
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sizes fear* inage analysis our unit areas represented a larg^ area on 
the ground than we would have likedi 

2 , Interpretation -problems 

It did not appear that we could satisfactorily disentangle the 
different and variable influences of picture processing and cloud 
Gharacteristics in order Idiat tje nii^t, with reasonable confidence, assess 
the cloud differences for their Statistical significance. Consequently 
we have expressed our interpretation of our results only in quali^tive 
terms. 


3. Schedule problem 

In view of tlie wide-ranging nature of our project we may be unable 
to devote as much time to the theme of this report as its theme deserves ^ 
However, we hope to apply the methods developed hane to seme further 
images for two purposes: 

a) To obtain further information on the relations between the 
urdt sizes of the areas measinred ty thn microdensitometer and . 
the brightness characteristics of the clouds and cloudy areas; and 

b) To test Hie similarity between such briglitness characteristics 
for similar- cloud fields vies^ed under contrasting solar illuraina.- 
tion conditions, e.g. around the summer and winter solstices. 


VI DATA QUAUry AND DELIVERY . 

Our inability to produce useful positives from the 70 mm negatives 
i^eived fron NASA continues to affect our vKrk adversely^ Test positives 
developed fron these negatives have been of a lower quality than Ihose 


we have produced by copying the 70 mm positives onto contact negatives 
and then developing these. We are not sure why this should be the case. 


Delivery of imagery secans to have slewed ftirtlrer as our project 







period has passed. We are still avEiting imagery for I'ferch 1976 at the 
tome of wfitiiig (mid-August 1976). This has had seme effects v^jon Ihe 
sequence of work in our study; some of these effects were alluded to in 
the Introduction to this Report. 


VII I^COMMEMDAFIQNS AND CONCXUSIOHS 

Although the assessment of the brightness characteristics of clouds 
imaged multispectrally by Landsat is a worthy topic of enquiry in its own 
right, much of our interest in this matter stems fron our broader work 
in the Applied Climatology Laboratory of the Department of Geography in 
the Urdyersity of Bristol. This is concerned with the search for an 
automatic (objective) method for cloud image analysis and cloud recognition 
in relation to a wide range of potential areas of application. The most 
diract methods of objective cloud mapping include: 

a) Brightness contour mapping, as practised vjith both visible 
ard infrared iiagety in the Applications Branch! of the 
National Oceanic and Atmospheric Administration, Canj) Springs, Md. 

b) Brightness/texture analysis, a basis for vSdch has been 


develcped in this DepaJTto (Harris ard Barrett, in preparation)} 
and c> Multispactral image analysis , summarised by Shenk, Holi±> and Neff 
■/ ^-(isTe), - . 


The first ot these would be difficult to develop into an operational 
cloud recognition method because of the ambiguif / wliich is evident in 
contour naps prepared frKmi either visible or infrared imagery Tshen 
relationships between’ image brightness and cloud tj,pe are considered. , , 

For example in visible imagery moderately thick low to middle level 
stratified clotd may be difficult to distinguish frem thick alto- or even 
cirrostratus ; these cloud types (and others) all reflect incident radiation 
strongly. In infrared imag( 2 y' deep convective cloud Can give hriglitness 
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teniperatures as low as those associated witii deep layered cloud. At bast 
such a method would be ari enhancement technique associated W3tn which 
thba?e woxild have to be sans qualitative assessment of the brij^tness 
contours in tcatis of cloud types. 

The second is much more attractive in that it enploys two criteria 
rather than one for the objective differentiation of cloud type. Ihitial 
tests suggest that brightness/ texture distributions for broad areas of 
cloud can be subdivided quite well into different cloud categories by 
discriiijiant analysis. The chief disadvantage of such a method arises 
from the coiitm sanpling iieth^ which subdivides the initial data array 
into sub-arrays (say 5 x 5 or 6 x 6 picture points apiece) in order that 
texture measures might be mads: texture is essentially a field characteristiG 
not a point characteristic as is the case with picture brightness . One 


si^iif leant consequence of tltLs is that seme of the resolution of the 
original data is sacrificed in brightness/texture analyses, on which 
Cloud reoDgnition. depends. Even then sane cloud types may be hard to 
separate, especially deep, id-despread stratiform/layered cloud and large 
convective dloud masses topped by dense anvils of cirrus. 

The third possibility v»uld appear tc have an advantage over the first 


in that a niinber of cloud responses would be available for each picture 
elenent, increasing the possibility that worthwhile cloud identification 
could be basod thereon. It would also appear to have sn advantage over the 
second in that the full spatial resolution of the original picture informa- 
tion obtained fixm the satellite could be used for cloud mapping and 


recognition. Hewever, the present Landsat MSS is not ideally suited for 



such ah q^xaration because of the narrowness of its coverage of the visible- 
infrared region of the elacrtxcniagnetic spectrum: the addition of a 


3 CA channel in the broad atmospheric window waveband in thermal infrared (as 

Is 


P-i planned for Landsat C) should significantly help to differentiate further 


SO 


the spectral signatures for different categories of cloud; the addition 
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of a furtJier channel deeper into the infrared region (c. 20~23 Mm) mi^rt 
help still rorare, especially in diffejrentiating between deep layered 
clouds and cumulonimbus topped by thick cirrus varieties. In their 
earlier work Barnes and Chang (1968), Greaves and Chang (1970) and Lo 
and Johnson (1971) concentrated upon rather lov; resolution weather 
satellite data when seeking tc asses the possibilities of multispectral 
cloud type identification procedures. Similarly Shenk and Holub (1973) 
and Shehk, Holub and Neff (1976> based their concliisions on MRIR analysis, 
and adndtted that their proposed Cloud Type Deci-sion Matrix might be 
ejqsected to t\) 0 !rk well without further assistance and/or refinement only 
over tropical oceans. We vrauld li]«:e to think that some future satellite 
syston might ecoploy a multispectral sensor system with a daily, global 
mapping capability and a resolution more in keeping with the present 
landsat MSS than the Tiros or Nimbus MEUiR - but designed for simultaneous 
Earth observation/atmospheric assessment purposes . It is probably only 
in this way that a multispectral cloud identification system could be 
Operated economically in the foreseeable future - econaity of effort as 


well as cost being significant when it is ramembered that a brightness/ 
texture cloud identification scheme is well vdthin the irecdiiis of practical 
possibility, this sufferir}g chiefly by cemparison on account of the 
resolution problem. 

It would seem that considerable benefits might accrue to all concerned 
were atmospheric and terrestrial inVesti^tors to collaborate more closely 
in future satellite design and operations; the problems of meteorologists, 
oceariDgraphers and Earth resource scientists are, perhaps, not cruy more 
closely interrelated then is ccnnonly acknowledged, but more amenable to 
solt±ions wMoh overlap considerably even though they might not precisely 


coincide. 
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